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Abstract. We discuss a recently proposed analytical formula for the eigenvalues 
of the Gaussian well and compare it with the analytical expression provided by the 
variational method with the simplest trial function. The latter yields considerably 
more accurate results than the former for the energies and critical parameters. 

In a recent paper Koksal [1] proposed a simple analytical expression for the 
eigenvalues of the attractive Gaussian potential 

V(r) = - 7 e- Ar2 (1) 

where 7 > is the well depth and A > determines its width. From the expansion of 
the potential-energy function about its minimum V(r) = —7(1 — Ar 2 + ^-r 4 — . . .) and 
perturbation theory the author derived an expansion for the energy of the form 

E nl = E»° + AE^ + AE%> + ... (2) 

where n = 0, 1, . . . and I = 0, 1, . . . are the radial and angular-momentum quantum 
numbers, respectively. The first term E^° is the sum of the minimum potential energy 
—7 plus the harmonic oscillation about this minimum. 

It is well known that the perturbation series (EJ) is suitable for sufficiently deep 
wells (sufficiently great 7) [2J. In order to obtain a better expression Koksal rewrote 
that perturbation expansion in terms of an exponential function. Since the author did 
not describe the general strategy clearly we conjecture that the main idea is embodied 
in the following expression 

E^ l = l(E^ + 1 )- 1 e-( E "°^^ (3) 

For large 7 we expand the exponential function and obtain the first term of the 
perturbation series (J2j) exactly and the approximation — ^- (^E^f + 7) to the second 
one. Koksal did not discuss the agreement between the analytical formula fl3]) and 
the perturbation series (j2J). Consequently, without further justification this expression 
can be considered to be an empirical formula and its validity determined solely by the 
accuracy of the results. It is worth noting that we can write equation ([3]) without 
recourse to perturbation theory because we only need the term of order zero. 
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Koksal carried out some calculations for the particular model parameters X = l/a 2 B 
and 7 = 400 Ryd, where as is the Bohr radius and Ryd the Rydberg energy. The 
approximate formula ([3]) appears to approach the numerical eigenvalues reasonably well 
for some values of the quantum numbers. However, we do not know the actual accuracy 
of the empirical formula ([3]) because the author did not report results for other well 
depths. What we already know is that the accuracy of the empirical formula decreases 
with I and most remarkably with n pp. 

The purpose of this comment is to test the accuracy of the empirical formula ([3]) 
more extensively and compare it with a simple analytical expression obtained by means 
of the variational method. 

The Schrodinger equation is 

Hip = E( 7 ,X)ip 

H = ~^ + v ^ < 4) 

where m is the mass of the particle which Koksal chose to be the electron. It is always 
convenient to work with a dimensionless eigenvalue equation that we easily derive in 
terms of the dimensionless coordinates r' = r/L, where L is an appropriate length 
unit. The Schrodinger equation thus becomes H'tp' = E'tp', where H' = mL 2 H/h 2 and 
E' = mL 2 E/h 2 . If, for example, we choose L = A -1 / 2 then we obtain 

H'= - ^V' 2 - £e~ r ' 2 

a (5) 

where V' 2 = L 2 V 2 . Note that E' depends only on the parameter £ since E'(£) = 
The dimensionless version of the empirical formula ([3]) is 

KY(0 = l(2n + i + 3 -)^- ^-K^+Dv^ (6) 

so that the discussion of its accuracy is greatly facilitated by the fact that we need to 
vary just one model parameter. Note that when A = l/a 2 B then £ = 2 2 "^b j s ^gjf 
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well depth in Rydberg units Ryd = 2 J^ 2 ■ Therefore, the particular values of the model 
parameters 7 and A chosen by Koksal correspond to £ = 200. 

In a recent pedagogical article Fernandez [3] discussed the application of the 
variational method to the one- dimensional Gaussian well (see also [1]). We can apply 
the same approach to the Gaussian well in three dimensions. Following those papers we 
choose the simple trial function 

v (r) = Nr l+1 e~ ar2 (7) 

where N is a normalization factor and a > is a variational parameter (we drop the 
primes on the dimensionless variables from now on). The optimal value of a is given by 
a root of d (H r ) /da = 0, where H r is the radial Hamiltonian 

Hr --2d^ + ^2^-^ e ® 
We thus obtain 
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We can proceed in two alternative ways: either, given £ we solve the first equation 
numerically for a and then obtain the energy or we obtain both £ and the energy 
analytically for a set of values of a (a parametric equation for the energy). 

Fig. [T] shows the eigenvalues E' ol for £ = 200 and several values of I calculated 
by means of equations fl6]) and (Q. The highly accurate eigenvalues provided by the 
Riccati-Pade method (RPM) [5] can be considered to be exact for present purposes. As 
discussed above £ = 200 corresponds to the potential parameters chosen by Koksal. We 
appreciate that E' f deviates from the exact result as / increases. On the other hand, 
the variational energy Eq^ deviates so less noticeably that it appears to agree exactly 
with the exact energy in the scale of the figure. 

Fig. |2] compares the approximate ground-state energies E' q , E' Q ar and the exact 
RPM ones for a range of values of £. We appreciate that the variational energy is closer 
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to the exact one for all values of £. However, the empirical expression appears 
to yield reasonable results for the ground-state energy for all those values of the well 
depth. Note that the largest potential parameter £ = 30 in Fig [2] is considerably smaller 
than the one chosen by Koksal. It is well known that the deepest the well the more 
accurate the results of perturbation theory [2] . For this reason the values of the potential 
parameters in Fig. [2] pose a good test for any formula based on perturbation theory. 

It is also well known that the Gaussian well supports a finite number of bound 
states and that there are critical values of the potential parameter £ for which bound 
states are exactly at the threshold of the continuum spectrum E = 0. In other words, 
there exists a bound-state eigenvalue E' nl provided that £ > Cnz, where E' nl (^ n i) = 0. 
Koksal did not discuss this important problem by means of the empirical formula ([3]) 
altough it is obvious that we can obtain estimates £^ from the roots of E' n f(£) = 0. 
We cannot solve this equation exactly but the numerical calculation is simple enough. 
On the other hand, from the variational energy (Q we obtain a^i = (2/ + l)/4 and the 
extremely simple analytical formula 



Fig. [3]shows ^qi, £gf r and the accurate numerical results obtained by Liverts and Barnea 
[6]. It is clear that while ^ merely follows the trend Qf 7 " is almost indistinguishable 
from the exact results in the scale of the figure. More precisely, the accuracy of 
decreases noticeably with I while Qf r remains remarkably accurate for all I values. 

Since the author did not give a clear justification for the empirical formula ([3]) nor a 
sound procedure that may be applied to other problems we assume that the sole purpose 
of the paper was to obtain an empirical formula for the eigenvalues of the Gaussian 
well. This assumption is supported by the fact that Koksal did not attempt to derive 
a similar expression for the eigenvalues of the Yukawa potential already treated by the 
same perturbation method in an earlier paper [7]. On the other hand, the variational 
method discussed above is not restricted to the Gaussian well and can be easily applied 



rvc 

so; 



var 



(2l + 3)~ 

21+ 

8(2/+ I) — 



(10) 



Energies for the Gaussian well 6 
to the dimensionless Schrodinger equation for the Yukawa potential 

V(r) = -^e~ r (11) 
By means of the trial function 

<p(r) = Nr l+1 e~ ar (12) 
we obtain the variational parametric formula for the energy 
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and the critical parameters are given by the simple analytical expression 

2 2l (I , ^2 1+3 

par _ L \ L I X J / u n 

Fig. H] shows the remarkable agreement between this formula and the accurate numerical 

results of Liverts and Barnea [B] • 

Finally, we summarize the main conclusions of this comment: 

First, Koksal's empirical formula is far less accurate than the analytical expression 

provided by the simplest variational function. It is true that Koksal's formula applies 

to states with n > while the variational method does not yield simple analytical 

expressions for such states (the Rayleigh-Ritz method suitable for them should be 

treated numerically). However, it is also true that Koksal's empirical formula becomes 

considerably less accurate as n increases pQ and here we have just compared the results 

for the most favourable case n — 0. 

Second, the variational method applies to other problems as we have just illustrated 

by means of the Yukawa potential. For unknown reasons Koksal did not attempt to 

apply the same approach to other models for which perturbation corrections are already 

available [2] as it is the case of the Yukawa potential [7]. 
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Figure 1. Energy eigenvalues E' Q f (red squares), E'^ ar (blue circles) and E' 
(solid line) for £ = 200 
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Figure 2. Ground-state energy E' q (red squares), E' () Q ar (blue circles) and E' 
(solid line) for a range of £ values 
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Figure 3. Critical parameters £q£ (red squares), Q^f r (blue circles) and £of° c * (solid 
line) 




Figure 4. Critical parameters Q® r (blue circles) and ^ act (solid line) 



